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One of the principal soiirccs of interchannel interference in mullichannel 
FM and PM systems is the dependence of the attenuation and phase shift 
of the transmission path on frequency. Here we study the interference pro- 
duced by a simple kind of such a dependence, namely that due to a single 
echo. Besides being important in themselves, the residfs are of interest be- 
cause they may be used to estimate the interference in other, more complicated 
cases. 

In this paper expressions are derived for the interchannel interference 
produced by echoes of relatively small amplitude. Several important special 
cases are studied in detail and curves that simplify the computation of the 
interference power are presented. 
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INTRODUCTION 

The systems we shall consider are of the FDM-FM and FDM-PM 
types; that is, systems in which the composite signal wave (the "base- 
band signal") from a group of carrier telephone channels in fieqiiency 
division multiplex (FDM) is transmitted by fref|iicncy modulation (FM) 
or phase modulation (PM). Such methods are currently being used in 
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the Bell System to send large groups of telephone channels by microwave 
radio. If the FM signal is a<x-ompanied by echoes, which may be due to 
reflections in the equipment or in the transmission medium, the wave of 
instantaneous frequency versus time is distorted in a nonlinear manner 
and interchannel interference occurs. Here we shall be concerned with 
this interference. 

The distortion has been analyzed in a number of previous publica- 
tions^"'' for the case in which the base-band signal may be represented 
by one or more sine waves. However, when the number of telephone 
channels is not .small, the sine wave representation becomes unwieldy 
because a large number of both low^ and high order modulation products 
must be considered. Here we avoid this difficulty, at the cost of some- 
what more complex analysis, by using a band of random noise to repre- 
sent the multiplex signal. 

It has been found in practice that such a random noise signal of appro- 
priate bandwidth and power adequately simulates a composite speech 
signal. For studies invoh'ing interchannel interference the energy cor- 
responding to some particular telephone channel is removed. When suclx 
a wave is impressed on the frequency modulator and the resulting FM 
wave is transmitted, detected, and fumlly demodulated, the received 
output in the originally clear channel represents interchannel mterfer- 
enee. Measurements of this type have been discussed by Albersheim 
and Schafer.^ 

One of the principal sources of interchannel interference is the vari- 
ation of the attenuation and phase shift of the transmission path with 
frequency. An echo produces a simple form of such a variation. In fact, 
it is often possible to estunate the effect produced by a more complex 
type of variation by comparing it with a roughly equivalent echo. 

It is the purpose of this paper to develop formulas for the interchannel 
interference produced by echoes of relatively small amplitude in FM 
and PM systems. General expressions are given in the form of definite 
integrals which may be evaluated by numerical integration. Approxima- 
tions are obtained in a number of cases of importance and numerical 
tables and curves are furnished to facilitate applications to specific prob- 
lems. 

We wish to express our thanks to Miss Mary Corr and Miss Barbara 
Fischer who have performed the rather lengthy computations required 
for this investigation. 

1. DEVELOPMENT OF GENERAL FORMULAS 

Since the tolerable amount of interchannel interference in multiplex 
telephony is small, the practically important situations are characterized 
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by echoes which are relatively small relative to the main transmitted 
wave. It is necessary to consider in detail only the case of one small echo. 
Effects of small multiple echoes may then be calculated by superposition 
of the effects of single echoes. Our distortion problem then reduces to 
the following: An original signal 

Eo(i) - Esm[pt + ip(0] (1.1) 

is received along with an echo 

E^iO = rEmn {p(t - T) + <p(t - T)] (1.2) 

Here r represents the ratio of echo amplitude to the amplitude of the 
principal received component, p is the unmodulated carrier freciuency, 
T is the delay difference of the paths and ip(t) is the phase variation 
caused by the multichannel signal. 

Calculation of the phase of the original signal plus its echo is equivalent 
to calculating the value of 6 in the representation 

sin X '+ '■ sin (x -\- y) = V sin (.i: + 6) (1 .3) 

with X, !/, r, V, and 6 real numbers. V and 6 are determined as functions 
of r and y by the two ecjuations obtained from the sin x and cos x por- 
tions of (1.3). When r « 1, V approaches unity, and 8 is proportional to 
r; direct expansion in powers of 6 shows that to a first order approxima- 
tion 

6 = r sin y (1.4) 

In our case, 

X = pt + <p{t) (1.5) 

y = ^(t - T) - <p(t) - pT (1.6) 

Hence the phase error produced by the echo is 

e(i) = rsm[v(f) - pT], (1.7) 

lit) = <p(l ~ T) ~ <p(t). (1.8) 

Our problem is to study the power spectrum of 6(t) when (p(l) is a random 
noise wave. 

The treatment will be restricted to the practically important case in 
which the noise source is free from dc and discrete sinusoidal components. 
A sufficient mathematical condition, which we shall adopt to insure this, 
is that the power spectrum of the noise have limited total fluctuation. 
We shall assume also that the noise wave is a member of a Gaussian 
ensemble and hence that all the statistical properties are derivable from 
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either its power spectrum or autocorrelation function. Which of the two 
quantities we use at any stage of the computations is a matter of relative 
expediency. We shall adopt a uniform notation in that an ensemble of 
time functions x{t) will be said to have a power spectrum Wx{f) and an 
autocorrelation function Rx{t), i.e., 

K.(r) = ave [xit)x(t + r)] = / w,if) cos 27rr/ d/ (1.9) 

W,{f) = 4 f RXt) COs2x/rdr (1.10) 

The power spectrum is proportional to the mean square of the response 
of a filter of bandwidth df centered at / when x(t) is applied as input. 

In the echo distortion problem, the power spectrum w^(f) is the quan- 
tity given and the power spectrum waif) is the one desired. The cor- 
responding autocorrelation functions R^{t) and Reir) are calculable from 
the corresponding power spectra and vice versa. The power spectrum is 
the more convenient choice as a function to compute when we deal with 
the noise response of a linear system with known transfer admittance 
Y(if), since the power spectrum of the response is merely [ Y(if) ]^ times 
that of the input. For example, differentiation is equivalent to transmis- 
sion through a transfer admittance i'lirf and hence multiplies the power 
spectrum by 4^^V^ The inverse process of integration divides the power 
spectrum by iir'Y- When a highly nonlinear operation is performed on the 
noise, or when the operation is a linear one more simply described in the 
time domain, the autocorrela,tion function may be simpler to compute. 

The first step in the solution of the problem is to evaluate the statis- 
tics of the noise wave ensemble, 

(.(0} = W(t- T) -<p{i)] (1.8) 

in tei-ms of the statistics of <p(0- This can readily be done in terms of 
power spectra by calculating the transfer admittance of a two-path trans- 
mission system .such as the right-hand member of (1.8) defines, or in 
terms of autocorrelation functions by averaging the appropriate time 
functions. The latter procedure gives : 

R,{t) = ave [vH) v(t + r)] 

= ave {{<p(t - T) -<p{t)]W{t - T + t) -<p(t + r)]} 

= ave [<p(t - T)ip{t- T + r)] + ave y{t) <p(t -j- t)] 

- ave y(t - T)<p(t - T -\- T -\- T)] 

- ave [tp(i) <p(t + T - T)] 

« 2fl/T) - R^(t + T) - R^(t - T) 
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In the averaging process we have made use of the facts that the average 
of the sum is equal to the sum of the a\'erages and that ip(l) is a stationary 
ensemble, allowing us to replace i — T by t when averaging. Since the 
two-path transmission system is linear and invariable, the i;(/) ensemble 
remains Gaussian. The corresponding power spectra relations are, from 
(1.10), 

«;„(/) = 4 M?//) sm' 7r/r (1.12) 

Our next step called for by (1.7) requires the evaluation of the auto- 
correlation function of a sine function of a band of noise. The solution 
may be obtained* from (3.2 — 7) of Reference 5, which gives a general 
theorem for a Gaussian noise ensemble expressible in our notation as 



ave |exp {iax{t) + ihx{t + t)]) 



1 , ,2 (1.13) 



= exp [ - ^^ RM - ohR^ir)] 

In this equation a and h are real constants. Multiplying both sides by 
the constant exp (iff), /3 real, and equating real parts gives the more 
directly applicable result: 

ave {cos [ax{t) + hx{t + r) + /3]1 

= exp [- -^^ fl.(0) - ahRM~\ cos H (^"^ 

Referring now to (1.7) we write 

i?s(r) = ave {r sin \v{t) - -pT] r sin [vit + r) - pT\\ 

I - 

J.2 . . ^ . 

= - ave {cos [v(t) ~ vit + t)]} ' ^ 

2 

- ^ ave {cos [vH) + u{t + r) - 2pT]} (1.15) 

-^exp[- /t'„(0)+ R.(t)\ 

- ^ exp [- RM ~ /?>(r)l cos 2pT 



* Tliis cipplicatinii of the characteristic functions method of attack has heen 
employed in similar situations by D. Middleton in Uefereiice and M. K. Zinn 
in unpublished memoranda. 
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Note that since the operation of taking the sine is nonlinear, the e{t) 
ensemble is not Gaussian and neither the power spectrum nor the auto- 
correlation is sufficient to give a complete statistical description. Either 
will be sufficient for our purposes however. 

The power spectmm of the distortion ensemble e(t) is now determin- 
able from the Fourier cosine transfoi-m of 4 Rb(t) as indicated by (1.10). 
The remaining steps are simplified, however, if we perform two pre- 
liminary operations on Reir) before the final Fourier transform is calcu- 
lated. 

We first observe that if Rb(t) contains a component C which does not 
vary with t, the ensemble d(t) contains a dc component 9 = \/C. The 
presence of C complicates the integration and hence we subtract out such 
a term before taking the Fourier transform. The value of C is given liy 

C = lim fie(T) 

r-»a) 

= i- e^'^"" lim [e""'^* - e"""^^' cos 2pT] (i.iq) 

= ^\-''.«>(l -cos2p3') 

In calculating the limit we have made use of the fact that, for our as- 
sumed i^oise wave, Rv(t) must approach zero as r becomes infinite. Our 
autoc((rrelation function of interest is thereby reduced to 

/2^_5(^) = Le-"'*"' {e"'''' - 1 - [e""'*^* - 1] cos2pr} (1.17) 

which is the autocorrelation function of the ensemble Q{1) - 9. 

Our second preliminary operation on Reir) is suggested by our ulti- 
mate goal, namely the calculation of the interchannel interference spec- 
trum i(v(/) (i.e., wjf) (If is the average power received in an idle channel 
of width df when the other channels arc carrying signals). We note that 
the spreading of the original spectrum into initially vacant frequency 
ranges occurs solely because of the nonlinear dependence of Re-eiT) on 
Rv(t). The disturbance received in an idle channel is therefore produced 
by this nonhnearity. The part of Rs-M which varies linearly with 
R„(t) may be expected to represent the linear transmission, and the 
difference between Ro-eir) and its linear portion to represent the non- 
linear transmission. In other words, subtracting the linear portion from 
R6--b(t) is equivalent to lemoving the linear transmission from any 
chaiinel without disturbing the nonlinear contributions from the re- 
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maiiiing channels. These considerations lead us to set the autocorrela- 
tion function fl,.(r) corresponding to wdf) equal to Rb--'b(t) minus its 
linear portion. The worlc of Appendix I shows that this e(|uality is 
rigorously true. 

In order to perform the subtraction it is convenient to write (1.17) in 
the form 

lie-eir) = F{R,{t)] (1.18) 

where in our ca-se for a general variable z, 

F{z) = ^e"''""'"k* - 1 - ie^' - 1) cos2p3'] (1.19) 

The portion of Ra-'air) which varies hnearly with Rv{t) is F'{0) Rv(t) 
where 

F'{z) ^ ^ e-'^^'^ie' + e~' cos 2pT) (1.20) 

F'm = '- c-^-'^'d + cos 2pr) (1.21) 

According to the foregoing discussion and Appendix I, the autocorrela- 
tion function and the power spectrum of the interchannel interference 

are given by 

Rc(r) - Rs^sir) - F'(0)RXr) 

- [e-"'^'^ - 1 + fi,(r)]eos2pT} 

.(/) - 4 /" /?,(r) cos2T/rrf7 (1.23) 

Jo 



w. 



The power spectrum of ^(0 — may be simply expressed in term.s of 
uv(/) without further integration. We find by applying the Fourier tran.s- 
form relationship of (1.9) and (1.10) to (1.22): 

we-eU) = 4 r \R,{t) -f- F'(0)7?„(t)] cos 2wfrdT 

JQ 

= Wcif) + F'iO)w,U) (1.24) 

= «'c(/) + 4f '(0)w',(/) sin' tt/T 

= WcU) + 2rV"^'"(l + cos 2pT)w^U) sin" irfT 
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where we have used (1.12) and (1.21). The power spectrum Wg(f) of 6 is 
obtained by adding to (1.24) a spire at/ = to represent the power in 
6, the dc portion of 9. 

The results obtained up to this point may be summarized as follows: 
Let ,p(t) be the phase variation produced by the impressed multichannel 
signal. The echo (1.2) produces a phase error d(t) (1.7) in the received 
signal. e(t) has a dc component 8 given by the square root of (1.16). 
9(t) — 9 is a function of tune which fluctuates about zero and has the 
power spectrum Wg_i(J) given by (1.24). Now consider the problem of com- 
puting the interchannel interference. One procedure would be to consider 
the case in which all but one of the channels are loaded. Then the power 
spectrum Wp(/) of ^(t) would have a narrow slot in it corresponding to the 
zero power in the unloaded channel. The values of v)e-i(f) computed from 
this w^(f) for values of / within the slot would give the channel inter- 
ference spectrum (for phase modulation). However, it turns out that as 
the slot width approaches zero, these values of We-ed) approach those of 
w,(f) where w,{J) is computed from (1.23) on the assumption that all 
channels are loaded. In other words the w^(f) used in the calculation of 
(1.23) has no slot. (Actually the same value of w^(f) is obtamed whether 
ivM) fias a slot or not. Of course, the slot must be vanishingly narrow). 

Almost all of the preceding work pertains to the power spcctriun of the 
phase error d(t). For the sake of completeness we shall give the power 
spectrum of the complete phase angle, Q(t) = (p(t) + d(t), of the output. 
In order to obtain all of the O(r^) term it is necessary to add another 
term to the approximation (1.7): 

dit) = r sin \v{l) - prl - ^ sin [2v{t) - 2pT] + O(r') (1.25) 

The autocorrelation function for Q(t) may be shown to be 

Ro(t) = R^(r) + ReCr) + ave [^(0 6(1 + r) + >p(t + t) 6(1)] 

= RM + Ur) - rR.(r) r-«"™"^ COS pT (1.2G) 

-\- r'R.(r) e-'"-''' cos 2pT + 0(/) 

The average has been evaluated by using 

ave \<p{t ± r) sin Wit - T) - M - vT]\ ^^ ^^^ 

= [RAT ± r) - i?^(r)l cos pT exp [R,iT) - K,(0)1 

and the corresponding result for 2ip and 2pT (Ri^r) = 4/?^(r)). Equa- 
tion (1.27) may be obtained from the result similar to (1.13) which con- 
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tains the three variables ip(t), tp{t ± t), ip{t — T) instead of x{t) and 
x{t + r). 

The power spectrum ivqij) of Q{t) may be obtained from (1.26) by 
replacing the autocorrelation functions in the coefficients by the cor- 
responding power spectra. By using (1.24) we may obtain the following 
expression for the power spectrum of the fluctuating portion Q{t) —Q of 
Q(0) Q =& being the dc portion of Q{i)'. 

Wq-qU) = loM) + wM) + ivM)[-r e-'^'''"' cos pT 

+ r e-""''' cos= pT + r' e-'"^''' cos 2pT] (1.28) 

The expression (1.12) for «Jd(/) shows that it is proportional to w^(f). 
Hence, assuming w^(/) to have a narrow slot corresponding to an idle 
channel, it is seen that (1.28) is in agreement with the fact that only wdf) 
contributes to the intcrchannel interference spectrum. 

We remark that for echoes produced by reflections in wave guides the 
angle 2pT in (1.22) and (1.24) usually contains many multiples of 27r. 
In most cases it would in fact be reasonable to average the effect of the 
term cos 2pT by assuming the angle to be uniformly distributed through- 
out 27r radians. This gives the result zero for the term since plus and 
minus values are symmetrically distributed. We shall, however, carry 
the term along in our formulas since it is of importance when the delays 
are small, as in multipath radio propagation, and when our results are 
used to estimate interchannel interference in general. 

2. APPLICATIONS TO FLAT NOISE SIGNAL 

In general Wp(/) depends on the type of preemphasis used in the chan- 
nel multiplex signal. Two representative conditions will be studied: (1) 
phase modulation (PM) in which channels at equal level arc impressed 
on a phase modulator, or the more usual equivalent situation of equal 
level channels which are differentiated before being impressed on a fre- 
quency modulator, and (2) frequency modulation (FM) in which chan- 
nels at equal level are impressed on a frequency modulator. The appro- 
priate power spectra are: 



PM: 


wM) = A , 


/.</</* 


(2.1) 


FM: 


wM) = n/4//. 


/.. <f<U 


(2.2) 



The latter form results because the phase is the time integral of the in- 
stantaneous fre{|uency, which has a flat spectrum. In PM, Po is expres.sed 
in (radians) '/eps, but in FM, Po is expressed in (radians/sec. )Vcps. In 
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the case of most common practical interest, the band of frequencies be- 
tween and /o is relatively narrow compared to the range /& — /„ . We 
have accordingly assumed that /„ approaches zero. In the FM case we 
cannot set /<, = immediately because the power spectrum would then 
become unbounded at the origin. It is found, however, that finite limits 
are approached for the actual quantities of interest which we compute. 
When the spectrum w^(f) of the signal has the form (2.1) or (2.2), 
expression (1.23) for the interchannel interference spectrum may be 
written as 

wM) = r\2irftr\G - II cos 2pT) (2.3) 

where wdf) df is measured in (radians) and 

G = 2e-"""' r [e""''' - R^r) - 1] cos audu (2.4) 

Jo 

H = 2e-^'™ r te"''"^^' + R.(r) - 1] cos au du (2.5) 

Jo 

u = 2irUr U = 2x^7 a = f/f, (2.6) 

ForPM: 

T^ f \ n r [2 sin w sin (u + U) sin (u - U)! ,„--, 

where Po/b is the mean square value, in (radians)^ of the PM signa^ 

For FM: 

E„(t) = A[-2(l - cos C/)cos u - 2uSi(u) 

+ (m + U)Si(u + f/) + (u - U)Si(u - U)] (2.8) 

In (2.8) Pofi is the mean square value (measured in (radians per second)") 
of the FM signal <p'(t). Consequently, A = (c/fbf where u is the rms 
frequency deviation, in cycles per second, of the signal. Equations (2.7) 
and (2.8) follow from (2.1), (2.2) and (1.11). 

When Rv{0) is small compared to unity, the same is true of Rv(t), and 
(2.4) and (2.5) lead to the approximations 

(-OS ' 

G ^ II w / Rv\t) cos au du (2.9) 

Jo 
The interchannel interference calculated from (2.9) is that produced by 
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the "second order modulation products" and is studied, together with 
other approximations, in Sections 3 and l. 

The quantity of interest in practice is the ratio of average interference 
power in the idle channel to average signal power in an adjacent channel. 
For PM the average interference and signal powers (in a narrow channel 
centered on frequency /) are the channel bandwidth times We(J) and 
wM, respectively. When the number of channels is large, the power 
spectrum does not change appreciably in going from one channel to the 
next and we may write 

^ = ^^ = 9-7^ f^ - ^^ '''' '^VT) (2.10) 

as the desired interference-to-signal power ratio. 

For FM the average interference and signal powers in a narrow channel 
arc (27r/y X (chaTuiel bandwidth) times wdf) and ^v^(f), respectively. 
These powers are measured in (radians per second)^ When we take the 
ratio Pj/Ps , the (2x/)^ X (channel bandwidth) cancels out and we have 
from (2.2) and (2.3) 



3. APPROXIMATIONS FOR G AND H PHASE MODULATION 

Table 3.1 contains various approximate expi'essions for the quantities 
G and // which enter expression (2.3) for the channel interference spec- 
trum Wc(J). The notation is explained in Section 2, and tlie results apply 
to the case in which w^(f) has the flat power spectrum (2.1). 

Case 1 gives the exact expressions developed in Section 2. Case 2 is 
the "second order modulation approximation", valid when RviO) « 1. 
Evaluation of the integral (2.9) for G with the help of 



£ 



sin u sin (u — U) , 

cos au au 



u{u - U) 






tU~' cos I ^ ) sin ( t/ - -^ ) , ^ a ^ 2 
0, 2 ^ n 



leads to tlie following expression for the quantify J appearing in Table 
3.1: 



612 
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Table 3.1 --Phase Modulation 



Case 
No. 


Restrictions on Parameters 
U aod Pifi, 


U and H 


Notes 


1 


No restrictions, 
wM) defined by 
(2.1) 


G defined by (2.4) 
H defined by (2.5) 
R,{t) defined by (2.7) 


Pofb = ave ynt)] 
Pi n 




P, 27rPoA 

{G - H cos 2p2'), 
U = 2wUT, 

a = f/h 


2 


2P.H (> - =V') 

«1 


G '^ H == 2w{P<,jbyJ, J defined 
by (3.1) 


"2nd Order Modu- 
lation" a jprox., 

J tabulate! iti 
Table 3.2, 10 log,o 
J plotted in Fig. 
5.1 


3 


U « 1 


G =^ H ^ 2T{P,hy 240 

.|12 - 30a + 200!= - a^] 


Special case of Case 
2 


4 


2Po/b« 1, 
I! » 1 


a = H ^ 2T(Po/t)= 

0-.!)0 + ^) 


Special case of Case 
2 

(? is a rapidly oscil- 
lating function of 
a 


5 


Po/6f/V3»l, 


G»H, 

G « (107r/P,)/,,f/=)''' 

exp (-10fiV4/'o/iX'^) 


When U « 1 
Case 3 applies when 

Pofb is small 
Case 5 applies when 

Pofb is large 


6 


t/»l 


f? = p-^o[/(bc , a) 

+ 2I(bi , a) cos all], 
H = c-*o|/(-6o , a) 

+ 2/{— 6i , a) cos «L'], 
fco = 2PoA , fei = -/'o/fi 
/(/), a) atiidiodin AppoiulixIII 


When Poh « 1. 
Case 6 reduces to 
Case 4. WhenPoA 
»1,G»//. Asfl 
increases, G oscil- 
lates between G+ 
and G^ given in 
Table 3.3. See 
Table 5.1 



J - 



a\/. , cos aU\ sin {2U - aU) 
-iK ^ 2 J^ 4f/ 



2 

--cos 



aU 

2' 



,_- 



(3.1) 



Values of J are tabulated in Table 3.2 for various values of a and k 
where U = /l:t/4 (or T = k/iSfb)). J is zero when a exceeds 2. 

Cases 3 and 4 in Table 3.1 follow directly from Case 2. In order to 
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Table 3.2- — Values of ./ 



k 


a = 


0.25 


0.50 


0.75 


1.00 


1.25 























1 


0.018 


0.009 


0.003 


0.001 


0.002 


0.004 


2 


0.227 


0.115 


0,041 


0.011 


0.023 


0.057 


3 


0.794 


0.434 


0.161 


0.049 


0,098 


0.231 


4 


1.500 


0.903 


0.352 


0.123 


0.250 


0.524 


5 


1.924 


1.2S4 


0.527 


0.224 


0.458 


0.810 


6 


1.5124 


1.385 


0.585 


0.332 


0.659 


0.936 


7 


1.712 


1.231 


0.505 


0.427 


0.773 


O.SOO 


8 


1.500 


0.994 


0.375 


0.506 


0.750 


0.494 


9 


1.M5 


O.SOO 


0.325 


0.579 


0.602 


0.228 


10 


1.245 


0.65S 


0.425 


0.657 


0.405 


O.ISO 


11 


1..107 


0.544 


0.643 


0.725 


0.262 


0.358 


12 


1.500 


0.472 


0.S83 


0.752 


0.250 


0.601 



obtain Case 5 we note that when f/ <K 1, expression (2.7) gives 

d~ sin u 



ESt) w -C/ToA 



du'- 



(3.2) 



The corresponding integrals for G and H could, if required, be investi- 
gated by the methods used to study Lewin's integral in Appendix III. 
However here we consider only the ease where U'Pofh is so large that (1) 

exp Rv(t) is tlie dominant term in the integral (2.4) for G, and (2) most 
of the contribution to the value of the integral comes from the region 
around ii = 0. The results for Case o then follow from (2.4) and the fact 
that (3.2) becomes 



R.ir) 



-'^"^(j-fo 



When U y> I, expression (2.7) shows that Rv(t) is small except when 
u is near or near U: 



Rv(- 



2 Pofb " ^ sin », 



/?„(t) ^ -P.A((( - O"' sin (u ~ U), 



u near 
II near U 



(3.3) 



These approximations, expression (2.4) for G, and the definition (A3-1) 
of I{h, a) give tiie results stated in Case 6. Ju.st as in Case 4, G is a 
rapidly oscillating function of a when U is large. It oscillates between 
G^ and G~ where 



G^ = e-''V(0,,a) ±2/(6, ,a)J 



(3.4) 

A'alues of G^ and G~ are given in Table 3.3 for various values of a = 
j/Ji, and Pafb [ill (radians)"^]. 
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Table 3.3 — Values of G^ and CT. The Upper Number of an Entry 
IS G^ AND THE Lower Number is CT. 



Pojb 
(radians)' 


a = 


0,25 


0.50 


0.73 


1.00 


1.25 


0.25 


0.384 
0.160 


0.338 
0.144 


0.292 
0.126 


0.245 
0.107 


0.197 

0.087 


0.148 
0.066 


0.50 


1.018 
0.504 


0.906 

0,464 


0.780 
0.415 


0.665 
0.357 


0.537 
0.291 


0.406 

0.222 


0.75 


1.55 
0.89 


1.30 
0.83 


1,22 
0,75 


1.04 
0,65 


0.845 
0.533 


0.645 
0.411 


1.00 


1.90 
1.22 


1.73 
1.15 


1.53 
1.05 


1.32 
0.92 


1.07 
0.76 


0.830 
0.596 


2.00 


2.16 

1.86 


2.04 
1.80 


1.88 
1,68 


1.68 
1.52 


1.43 

1.31 


1.17 

1.07 


4,00 


1.59 
1.57 


1.56 

1.54 


1.50 

1.48 


1.40 
1.40 


1.28 
1.28 


1.14 

1.14 



4. APPROXIMATIONS FOR G AND H — FREQUENCY MODULATION 

The various (.'use-s which we shall coiiyider for FM are roughly similar 
to tho.se cousiclered iu Section 3 for PM, and are listed in Table 4.1. The 
power spectrum w^{f) is assumed to be that given by (2.2). As pointed 
out in Section 2, the average FM signal power in a narrow channel of 
width A/ centered on frecineney/ is as.sumed to be PoAf (radians/second)^ 
if < / < h and zero if /i < /. The average interchannel interference 
power is (2x/)" Wr(/)A/ (radians/second)'. When / > U this gives all of 
the power present in the frequency interval A/. 

Case 1, Table 4.1, gives the exact expressions for G and //, and Case 2 
(corresponds to the "second order modulation" approximation which 
holds when K„(0) « 1, R,-{^) being computed from (2.8). iC is a function 
of a = f/Jb and U = 'IirfiT w^hich may be obtained by writing the in- 
tegral (2.9) as 

G ^ irh r R'(r)e"''' dr = 4-VA r idXx)w.,U - ■^) dx (4.1) 

where Rv(t) and Wv(J) are taken to be oven functions. From the defini- 
tions (1.12) and (2.2) for w^f) and w^(f) we have 



w„{f) = 



\PMr' sin' TfT, 
[ 0, 



I/I <.h 

l/l>A 



(^.2) 
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Table 4.1 — Frequency Modulation 


Case 

No. 


Restrictions on 
Parameters U and A 


Cand n 


Notes 


1 


No restriction, 
wM) defined by 
(2.2) 


G defined by (2.4) 
H defined by (2.5) 
R„ir) defined by (2.8) 


.4 = (Wft)\ 
IT ~ nns 
tdp(()/d(l/2ff, 
(J- = rms frequency 
deviation of 
signal meas- 
ured in cy- 
cles/sec. 
Pi /Pa 

= [rV/(2irA)] 
■[G - H cos 2pT] 


2 


2A[VSiiU) - 1 
4- cos U] « 1 


G ^ H '^ 2irA^a-WK 
K defined by (4.5) and (4.6) 


"2nd Order Modu- 
lation" approx. 

K tabulated in 
Table 4.2. See 
Fig. 5.2 


3 


AU^ « 1, 
U «l. 


G ^ H ^ 7rA2(7n2 - a)/i, 
£ o g 2 


Special case of Case 

2 
UK « a"t7H2 - (i)/8 


4 


A(7t « 1, 
f/»l 


G ^ H = 


2n■MH/a-^ < a < 1 
ttU^C/, a = 1 
0, fl > 1 


Special case of Case 
2 


5 


U <Jcl,Lewin'scase 


G =« e-Hib, a) 
H = e-'Ii-b, a) 
h = AU-' 


Case 5 agrees with 
Case 3 when b 
« 1. G » H for 
6 » 1. I{b, a.) de- 
fined and tubu- 
lated in Appen- 
dix III. See Fig. 
5.4. 


6 


4A » 1, 
U»l 


G ^ kf/,/(A sinlM/,)!"' 

exp [-2A (cosh y^ — 1)] 

j/i defined by 

>i r"' sinh 1! , 
2^1 = 1 " '" 


Value of G is in- 
dependent of f/. 
See Fig. 5.3. 


7 


U ^ 1, 
B = A[l-U~'smV] 


G = [;r/B]"^exp(-«V(4B)] 


(?»H 
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where the lower limit /a of the frequency band is taken to be very close 
to zero. When / > 2ft, the value of (4.1) is zero. When we take < / < 
2fb the limits of integration in (4.1) are x = f — fb and x = ft, ■ Changing 
the variable of integration in (4.1) from x to y = 2TrxT converts (4.1) 
into 

^icAV £^ y-\a - y)-' sm' | sin= {^) dy (4.3) 

where a = 2-irfT = aC/. By partial fractions 

aY\a - y)-' = y~' + ^a'Y' + (« - ?/)"' + 2a-\a - y)"' (4.4) 

Considerations of symmetry show that the a ~ y terms on the right 
contribute the same amount to (4.3) as do the ij terms. When the y~^ 
term is converted into y~^ by an integration by parts, (4.3) may be ex- 
pressed in terms of &i{x) and Ci{x) functions. In this way it may be 
shown that the approximation (4.3) for G and H has the value 

2-KA^V^KIa- = 2ir A'UK/a^ 

where 

K = {-U~' + /3"')(1 - cos C/)(l - cos^) 

+ {Si U - Si /3)(1 + cos a - 2a"^ sin a) 

+ {Si 2U - Si 2/3)(- cos a + a~' sin a) (4.5) 

-{- {Ci2U - Ci2\^\ - CiU -\- Ci\0 |)(sm a + a~' cos a) 

+ a-\2 + cos a)[loge {U/\ \) - Ci U + Ci \ \] 

Here a = all and /3 = a — C/ = (a — 1)[/. When / = 2fb,a has the 
value 2 and K is zero, as it should be. When f = fb , i.e., when a = I 
and a = U, 

K = (1 + cos U ~ 2tr' sin U) Si U 

+ (- cos f/ + U"' sin U) Si 2U 

+ {Ci 2U - CiU - loge 2) (sin U + U~' cos U) 

+ U-' (2 + cos t/)(iog« U + .577.. - Ci U) 

Values of K arc tabulated in Table 4.2 for various frequencies and de- 
lays (a = j/Jb and h = -HJ/t = 8/^7'). Cases 3 and 4 in Table 4.1 show 
that when U is very small K ^ a'U^{2 — a)/8 and when U is very large 
K ^ TT, 7r/2, or according to whether a<l,a=l,ora>l. 

Case 3 in Table 4.1 is a special case of Case 2 which may be obtained 
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Table 4.2 — V 


ALUES OF K 




k 


„ =0 


0.25 


ii.?n 


O-IS 


l.'KJ 


1.25 








0.0 


0.0 


0.0 


0.0 


0.0 


1 


<) 


0.006 


0.022 


0.041 


0.05S 


0.068 


2 





0.048 


0.164 


0.305 


0.422 


0-473 


3 





0.141 


0.490 


0.890 


1.20 


1.20 


4 





0.286 


0.961 


1.74 


2.19 


1.45 


5 





0.453 


1.57 


2.64 


3.05 


2.62 


6 





0.673 


2.16 


3.36 


3.41 


2.45 


7 





0.897 


2.69 


3.68 


3.13 


1.73 


8 





1.14 


3.13 


3.65 


2.40 


0.883 


9 





1.39 


3.45 


3.34 


1.62 


0.370 


10 





1.66 


3.68 


2.95 


1.20 


0.378 


11 





1.93 


3.79 


2.64 


1.32 


0.741 


12 





2.20 


3.79 


2.51 


1.89 


1.10 



by letting U become very small in (4.3). The value of Pi/Ps correspond- 
ing to Case 3 has lieen given by Albersheim and Schafer.^ Case 4 may be 
obtained by letting U become large in (4.5) and (4.6). 
When U is very small, expression (2.8) for Rv(t) becomes 



R^(t) ^ AUSi ^ sin u 



(4.7) 



Assuming yl C/' « 1 and substituting (4.7) in the "second order modula- 
tion" approximation (2.9) for G and H gives us another dei'ivation of 
Case 3 (sec (A3-2)). However, if the rms frequency deviation of the signal 
is so large that .4. U^ is not small, even though U « 1 , ^\-e have Case 5, 
the case investigated by Lewin." The formulas given in Table 4.1 arc ob- 
tained when (4.7) is set in the integrals (2.4) and (2.5) for G and H, and 
the results compared with the definition (A3-1) of 7(6, «). 
When U is very large, expression (2.8) for Rv(t) becomes 

[AWU - 2 cos u - 2uSi(u)], Q < n < U 

R.(r) ^ (4.8) 

[ 0, U <n 

Substituting (4.8) in (2.9) and integrating by parts t^-iee leads to another 
derivation of Case 4. The expression for G given in Case 6 is obtained 
from (4.8) and (2.4) by the method outlined in Appendix II. AwU is 
assumed to be so large that most of the contribution to the value of the 
integral (2.4) for G comes from the region around it = 0. It is also as- 
sumed that Rv(t) + 1 is negligible in comparison with exp Rv(t) in this 
region. This leads to the approximation 



LI F^ ^ / t 

Jo 



cos au (hi 



(4.9) 



which holds when t/ » 1 and AttU » 1. 
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The expression for G in Case 6 is merely the leading term in the asymp- 
totic expansion arising from the saddle point at u = iiji . When further 
terms in this expansion are obtained (using, for example, e(iuation 
(10.4) of Reference 7) it is found that the expression for Case (i should be 
multiplied by 

I [-Ics - oiji + s-(yi~^ - 2yi)] , „ 

"^ 48.4s^ -r •-■■ K • J 

where c and s denote cosh (/i and sinh i/i , respectively. The next term 
consists of 1/.4* times a function of //i , and so on. When iji becomes 
small, as it does when A becomes large or a becomes small, (4.10) becomes 

^ 48A 13824-4^ 

However, comparison of Tables 4.3 and 4.4 shows that the formula of 
Case G gives fairly reliable values of G when -4 is as small as 0.5 (f' must 
be large, of course). 

When (1 - a) is small and .4 « 1, but U still large enough to make 
.l7rf/» 1,(4.9) gives 

T + 2 arc tan ^^^^1 + QiA'/a) (4.11) 

At _\ 

This may be obtained by letting A become small in 

G^4iA/aY r e'"'Si{u)Fdu (4.12) 

Jo 

?/ = 1 — COS u — uSi(u) 

F - -2 cos auSiiu) + Si[{l - a)u] + iV[(l + a)u] 

whii-h may be obtained from (4.9) by integrating by parts twice. 

The formula for G given in Case 7, Table 4.1, is obtained when L' is 
taken to be of order unity and A is assumed to be so large that only the 
exponential term in the integrand of (2.4) is of importance. Most of the 
contribution comes from around u = where 

R,.(r) = R,{0) - ((14(1 - U~' sin [/) + • • • . 

When U» I, and .4 » 1, Cases 6 and 7 both give 

(7^(xA4)'"'exp[-aV(4A)l 

which leads to an e.xpression for Pi/Pa similar to one given by Alber- 



G ^ Aa 
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Table 4.3 — Values of G and H Obtained by Numerical 
Integration 





a = () 


(I.2.S 


0.50 


0.75 


Mm 


1,25 


a for r = 


3 














A = 


0.125 


0.6iW 


0.651 


0.590 


0.466 


0.326 


0.224 




0.25 


1.53 


1.46 


1.29 


1.06 


0.808 


0.561 




0.50 


2.17 


2.09 


1.97 


1.63 


1.38 


0.993 


//for r = 


3 














.1 - 


0.125 


0.424 


0.393 


0.335 


0.260 


0.182 


0.111 




0.25 


0.574 


0.52S 


0.447 


0.339 


0.233 


0.136 




0.50 


0.283 


0.255 


0.211 


0.156 


0.104 


0.039 


a for (' = 


6 














.1 = 


0.125 


2.62 


2.42 


1.S5 


1.18 


0.631 


0,257 




0.25 


3.26 


3.05 


2.50 


1.81 


1.17 


0.695 




0.50 


2.55 


2.46 


2.21 


1.86 


1.47 


1,10 


H fur U = 


6 














.1 = 


0.125 


0.802 


0.705 


0.465 


0.227 


0.067 


0.0031 




0.25 


0.266 


0.231 


0.148 


0.064 


0.0126 


0,0035 




0.50 


0.0092 


0.0079 


0.0048 


0.0018 


0.00017 


0.00019 




A = 0.0625 


0.1 2S 


0.25 


o,.>n 






a for (( = 


1 














f' = 


1.5 


0.0126 


0.045 


0.155 


0.446 








3 


0.111 


0.326 


0.808 


1.38 








a 


0.245 


(1.631 


1.17 


1.47 








12 


0.300 


0.657 


1.16 


1.47 







sheim and Schafer^ for long delay. When U « 1 and AU^ » 1, Cases 5 
and 7 both give 

G ^ {QtA-' ir-y" exp [-3a' .4"' U'-/2] 

Before the approximationH listed in Table 4.1 were developed, a number 
of value.s of G and H were obtained from (2.4) and (2.o) by numerical 
integration. These values, given in Table 4.3, are the best we have and 
may be used to cheek the various approximations. 

As an example of the values given by our approximations we take 
the ease U = G. In Table 4.4, "G — Case 6" has been computed from 
the formula gi\'en in Case fi, Table 4.1 (which assumes f — » %). When 
these values are compared with the corresponding ones In Table 4,3, it 
is seen that the agreement is not good for .4 = 0.125. Better agreement 
is shown by "G — Improved Case b" in which the \-aIues are computed 
from (2.4) and (4.8) by the method of Appendix II. It ditfers from "Case 
fi" in that F(u) of (A2-1) is exp [^^(t)] - R„(t) - 1 instead of merely 
exp [RM)]. 
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Table 4.4 — Approximate Values of 


G FOR 


[/ - fi 






fl = 


U-2,i 


u.Sd 


(),7.i 


1.00 


i.2.i 


G— Case 6 

A = 0.125 
0.25 
0.50 

G — Improved Case 6 
A - 0.125 
0.25 
0.50 


5.01 

3.54 
2.51 

2.62 
3.17 
2.50 


4.13 

3.26 
2.42 

2.37 
2.94 
2.41 


2.54 
2.58 
2.17 

1.74 
2.39 
2.17 


1,32 

1.80 
1.82 

1.06 
1.72 
1.82 


0.63 
1.14 
1.43 

0.56 
1.10 
1.42 


0.28 
0.66 
1.06 

0.26 
0.66 
1.06 



5. INTERCHANNEL INTERFERENCE POWER , 

The values of Pj/Ps , the ratio of the interchaiinel mterference power 
to the signal power, may be computed from the formulas (2.10) and (2.11) 
when G and H are known. One ^vould lii<e to have curves giving Pi/Ps 
for representative combinations of echo delay, signal power, and channel 
position which are likely to occur in practice. Howe\'er, the large number 
of such comlnnations coupled with the difficulty of computing G and H 
leads us to restrict ourselves mostly to curves for Cases 2 and 6 ia Tables 
3.1 and 4.1. In all cases the signal power Ps (percps) is taken to be equal 
to the constant value Po (measured in (i-adians)Vcps for PM and in 
(radians/sec)Vf'ps for FM) over the signal band (0, /&), and is zero out- 
side this band. 

Case 2 is the "second order modulation" approximation which, roughly 
speaking, applies when the echo delay is very short or when the rms devi- 
ation of the phase angle (for PM) or of the frequency (for FM) is small. 
Case 6 applies when the echo delay is very long. 

(a) "Second order modulation" approximation for PM — Table 3.1, 
Case 2. Since G ^ H, equation (2.10) may be written as 

10 logio (Pi/Ps) ^p + D,+ Ih + D, 

p ^ 10 login r~ 

Di = lOlogio (1 - cos2pr) (5.1) 

Z)2 = 10 logir, (P^fb) 

A - 10 logio J - - 

where p is the reflection coefficient expressed in decibels, Dy is a quantity 
which varies rapidly with T and whose representative value is zero. J is 
the function of a and U defined by equation (3.1). The quantities Pq , 
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/[, , o, U are defined by equation (2.1) and Case 1 of Table 3.1. Pafh is 
the average signal power in (radians)", a = j/Jh gives the channel posi- 
tion and U = 2TrfbT measures the echo delay. 

The approximation (5.1) holds when 2/*o/6(l — C/~' sin U) is small in 
comparison Axith unity. 

Fig. 5.1 shows Z>3 plotted as a function oi JbT for various vahies of a. 
The values of J which were used were taken from Table 3.2. Values of J 
for C/ « 1 may be obtained from the expression given for G in Case 3, 
Table 3.1. Case 4 gives anotlier special case. 

{h) Large delay, U» 1 for PM — Table 3.1, Case 6. When V is very 
large, G is a rapidly oscillating function of a. When in addition P^)fh is 
small, Case 4 shows that J fluctuates between (1 — a/2)/2 and 
3(1 — a/2)/2. The corresponding fluctuations in D2 are noticeable in 
Fig. 5.1 for the larger values of /iT. If 2Po.f6 is large compared to unity 
(and the delay is large), equation (5.1) no longer holds. In this case G ^ 
H and we may write (2.10) as 

10 logio (P./Pa) ^p + D, 

p = 10 logio / 

D, = 10 logio (2TPoftr'G (5.2) 

G ^ e~^V(lfo . ft) + 2I(b, , a) cos aU] 

where bo and 61, are defined in Case 0, Table 3.1. It is seen that as a in- 
creases from to 1, Di cscillatos rapidly between limits Di^ and Di' cor- 
responding to G^ and G^ which are defined by equation (3.4) and tabu- 
lated in Table 3.3. Table 5.1 gives values of Di'^ and D^" computed from 
Table 3.3. 

The entries corresponding to the values 0.25 and 0.50 for P(Jb must 
be used with caution in equation (5.2) since they do not satisfy 2/*o/b » 1 
and H is not negligible in comparison with G. 

(c) "Second order modulation" a-p-proximaiion for FM — Table 4.1, 
Case 2. By making use of the expressions for G and H given in Case 2 we 
may write equation (2.11) as 

10 logio(P//Ps) ;^ p + Z)i + Do' + D,' 

D^ = 10 logiD A (5.3) 

T){ = 10 logto UK 

Avhere p and Di are defined by (3.1) and .4 by (2.S). K is the function of 
a and U defined by (4.5). A is proportional to the signal power: A = 
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Fig. 5.1 — For "Second Order ModuUit.i(Hi" the ratio I'l/Ps for PM depends 
upon Ds as shown by equation (5.1). 



Table 5.1 — Values of Da^ and Da . The Upper Number of an 
Entry is D^ and the Lower One is Di' 





a = 


0.25 


0.50 


0.75 


1.00- 


PoU 


(radians)' 
0.25 


-6.12 
-9.92 


-6.68 
-10.37 


-7.31 
-10.96 


-8.07 
-11.67 


-9.01 
-12.56 




0.50 


-4. 89 
-7.94 


-5.40 
-8.31 


-6.00 

-8.79 


-6.75 
-9.45 


-7.67 
-10.34 




0.76 


-4.83 
-7.24 


-5.31 

-7.56 


-5,87 
-8.01 


-6.57 
-8.63 


-7.47 
-9.46 




1.00 


-5.20 
-7.12 


-5.61 
-7.38 


-6.13 

-7,77 


-6.78 
-8.33 


-7.67 
-9.17 




2.00 


-7,65 
-8.30 


-7.89 
-8.44 


-8.25 
-8.74 


-8.74 
-9.17 


-9.44 
-9.82 




4.00 


-11.99 
-12.05 


-12.08 
-12.13 


-12.25 
-12.31 


-12.55 
-12.55 


-12.94 
-12.94 
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{<^/fbf where o- is the rms frequency deviation of the signal in cycles per 
second. !)■< and Di play similar roles in (5.3) and {5.1). 

The approximation (5.3) holds when' 2A[U Si{JJ) — 1 + cos U\ is 
small in comparison with unity. 

The values of K given in Table 4.2 lead to the curves for D^' shown in 
Fig. 5.2. When I' « 1, Case 3 shows that D-,' ^ 10 logiu [a'r\2 - a)/8], 
and Case 4 shows that when f >> 1 (provided a < I and AwV « 1) 
D3' W 10 log,,, (tC) - 12.95 + 10 logio/ft?'. 

(rf) Large delay, U » 1 for FM — Table 4.I, Case 6. It has just been 
pointed out that when U becomes very large, P1/P3 depends upon the 
delay only through the term D3 ^ 10 logm irU (neglecting the rapidly 
varying term Di) if .1 fV « 1. If AUw » 1, Pi/Ps becomes independent 
of C/ as (/ ^ X . This follows from the fact that the formulas of Case 6 
allow us to write (2.11) as 



10 login (Pj/Ps) = p + D,' 
Di ^ 10 logio [Ga'(2TA)~'] 



(5.4) 
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Fig. 5.2 — For "SpcoihI Order Modulation" the ratio P j/Fs for FM dppeiids 
on Di' as shown by equation (5.3). 
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where p is defined by (5.1) and G depends only on a and A through 

G w Wyi/A sinh y^'"^ exp [—2A (cosh yi - 1)1 

_a_^ f"' sinh V ^^ ^5-S) 

2A h V ^ 

Fig. 5.3 shows D/ plotted as a function of A for various values of a. 
It is assumed that ArU » 1. 

(e) Sviall delay and large nns frequency deviation jor FM — Table 
4.1, Cases 5 and 7. It turns out that Case 5 (Lewin's, case, U « 1) and 
Case 7 (^ » 1 and U of order unity) may iae combined into a single case 
by taking the quantity h in the formulas of Case 5 to be 6A(1 — IT 
sin U) instead of AU'^. When V « 1, Case 5 is obtained. When A » 1 
the asymptotic expansion for 7(b, a) leads to Case 7 if f/ is 0(1). 

In order to put this combined case in a form suited to calculation we 
write (2.11) as 



Ft rVG(, H „ J\ 

~ = T- ( 1 - 77 COS 2pT I 

Ps 2tA \ G "^ } 



r^6(l-^-^ sin L/)aV^J(b.a) { U 

2-Kb \ G ^ 



(5.6) 



(5.7) 



10 logio {PilPs) f^ p + Di' + Ds' + A' 
where p is given by (5.1) and 

D{ = 10 logio (1 - (H/G) cos 2pT) 
D,' = 10 login (1 - U-' sin U) 
D,' - 10 logio 6a' e"' I{b, a)/(2rh) 
6 = 6A(1 - U~' sin U) 

Fig. 5.4 shows values of De, computed from the values of /(&, a) given 
in Appendix III, plotted as a function of b for various values of a. The 
maximum value of 3 db for Di occurs when AU~ <K 1 and cos 2pT = 
— 1 . When A U^ is large Di is approximately zero. 

Fig. 5.5 and 5.0 show, in a rough way, the regions in which the various 
approximations apply. For PM, the delay and the rms phase deviation 
(measured by C = 2x/b2' and (Pofbf', respectively) are the parameters 
which determine the type of approximation to be used. The regions in 
the [(Pp/fc),^'^ U] plane shown in Fig. 5.5 are marked with the numbers 
2a, 3, 4, 5, 6b where the integer indicates the case number in Table 3.1 and 
the letters a and b refer to Cases (a) and (b) in this section. Fig. 5.6 is 
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Fig. 5.3 — For long delayed echoes the ratio Pj/Ps for FM depends upon 
Di' as shown by equation (5.4). 



the corresponding figure for FM. The coordinates are U and A^'', where 
A^'^( = a/fh) measures the rms frequency deviation. The region numbers 
are 2c, 3, 4, 5, 6d, 7e where the integer.? indicate the rase number in Table 
-i.l. It will be noted that there are regions where no approximation is 
available. However, an answer may always be obtained by numerical 
integration of equations (2.4) and (2.5) for G and H. 

Fig. 5.7 shows values of Pi/Ps for the top channel (0=1) where the 
interference is often at a maximum in an FM system. The coordinates 
(A''^, fhT) arc essentially the same as those of Fig. 5.6. In order to 
simplify the plotting, the phase angle 2pT is assumed to be such that 
cos 2pT is zero so that the contours are given by 

The contours have been obtained in part from the various approxima- 
tions where applicable and in part from values obtained by numerical 
computation from the exact expression, While there are, of necessity, 
some areas of uncerl^ainty in Fig. 5.7, it should be adequate for most 
engineering purposes. No corresponding curves have been computed for 
the case of phase modulation. 
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delay the ratio P1/P3 for FM depends upon Dt' as shown by equation (6.6). 
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Fig. 5.7 — Contour.^ of constant interference in the top channel of a multi- 
channel FM system. 
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It should be noticed that Fig. 5.7 is plotted for the case a = f/fi, = 1. 
If Fig. 5.7 were plotted for values of a slightly less than unity there would 
not be much change except in the upper left hand corner (A small and 
JbT large), where the interference would tend to be 3 db stronger. This 
discontinuous behavior as a passes through unity is shown by Case 4, 
Table 4.1 where U » 1 and AttU « 1. When U » 1 and AtU » 1, as 
occurs when fbT~^ x> (with A held fixed), equation (4.11) gives 

10 log,o ^^^ « 10 log. [I + I arc tan (l^)] 

when ^ « 1. This shows that at /bT = f», the discontinuity arises for 
values of A near (1 - a)/ir. When A » 1 and U » 1, we have 

Pi .1 . .r.^ 2Tr"' 4.34a' 
10 logic ^^ f^ -15 logio A - 10 logio —^ jj- 

which changes only slowly as a -^ 1 . 

G. USE OF EQUIVALENT ECHO TO ESTIMATE INTERCHANNEL INTERFER- 
ENCE 

When a steady sinusoid exp (icot) is applied to a transmission medium 
of the sort we have under consideration, the output is exp (iut — a — 
1/3) where a and /? are the attenuation and phase shift, respectively. 
Distortionless transmission occurs when a is constant and )3 has a con- 
stant slope over the essential range of frequencies. Departures from these 
ideal c6nditions cause interchannel interference in multichannel FM and 
PM systems. Our evaluation of the mterference caused by a small echo 
may equally well be regarded as an evaluation of interference for a par- 
ticular kind of amplitude and phase distortion, namely that given by 

« ^ -r cos wT, ^^r sin wT. (6.1) 

These expressions are obtained by writing exp (ioil) + r exp [iu}{t — T)] 
in the form exp {-iojt — a — 7/3) when | ?■ | « 1. 

The analysis given by E. D. Sunde in Section 1 of Reference 8 shows 
that a minimum phase system in which /3 = r sin oiT also has a = 
— r cos 0)7' as in (6.1). This suggests a procedure for calculating inter- 
channel interference from phase data alone when (1) the distortion is 
known to be of the minimum phase type and (2) the variation of phase 
with frequency can be approximated by a sine function. In such cases 
we can apply our echo analysis directly by identifying r as the amplitude 
of the phase oscillation and T as the reciprocal of the period. 
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In carrier multiplex systems the sinusoidal approximation need hold 
only in the region around the carrier frequency /o where 27r/o = p and p 
is the radian frequency appearing in equation (1.1). For FM we shall, 
in this section, arbitrarily take the region to extend from/o — 4ff to/o + 
4ff where o- is the rma frequency deviation of the signal. For PM and the 
signal power spectrum given by equation (2.1) we may take the region 

A special case occurs when the nonlinear portion of )3 may be repre- 
sented as 02(/ — /o)V2 in the region of interest. We can think of r sin wT 
as going through several oscillations between / = and f = fo , and 
that a maximum, if a2 < 0, (or a minimum , if aa > 0) of r sin (^T occurs 
^t / = /o J ie. at 0} = p. The band of interest is taken to be narrow 
enough to lie in the immediate vicinity of the maximum. This sort of 
curve fitting is permissible since equation (2.3) shows that the inter- 
channel interference depends on the carrier frequency only through the 
term cos 2pT. Furthermore, constant terms and terms linear in frequency 
in the expression for /3 do not affect the amount of interchannel inter- 
ference. 

At the maximum mentioned in the preceding paragraph w = p, 
sin pT = I and cos 2pT = — 1. Near this maximum r sin coT is 

r cos 2t(/ ~fo)T^r- r4ir'(f - /o)'!rV2 (6.2) 

In order that this approximation may hold over the region /o ± 4a- (for 
FM), we require 

27r(4fr)r S 1 

We take T to be as large as possible, namely 

T = I/Stto- (6.3) 

in order to make r as small as possible, since our work assumes r « 1. 
Comparison of a2(/ — /o)V2 y^ith (6.2) gives 

r = -02(2x7")"' = -IQa^tx^ (6.4) 

which must be small compared to unity if our results are to be used. This 
result holds for 02 < 0. When 02 > expression (6.3) still holds for T 
but now r = \Qa2a'. Similar expressions hold for PM. These values of r 
and T may now be inserted in our formulas to determine the interchannel 
interference. 

From the definitions of A and U it may be shown that (6.3) is equiva- 
lent to AU^ = yiQ. Therefore the "second order modulation" approxima- 
tion given by Case 2 in Table 4.1 may be used. Also G — H cos 2pT is 
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approximately equal to 2G because cos 2pT = - 1. It turns out that the 
second order modulation approximation may also be used in the PM case. 
When a> is sufficiently small, considerations such as those above show- 
that the ratio of the interference power for j3 = a^if — Uf/2 (radians) 
to the received signal power at the frequency / is 

Pr/Ps = (a2amf a\2 - a) for FM 

P,/Ps = {a,f,y2f{P,U){\2 - ma + 20a' - a«)/30 for PM 

Here a = f/fb where (0, /6) is the frequency band of the signal. The rms 
frequency deviation of the signal for FM is <t cps, and the rms phase 
deviation of the signal for PM is {PofbY'' radians. The first equation in 
(6.5) comes from a special case of Case 2, namely, Case 3 of Table 4.1, 
and requires the additional assumption /f,/4o- « 1 (corresponding to 
(7 « 1). The second equation reciuires a similar additional assumption. 

Appendix I 

DERIVATION OF A GENERAL THEOREM ON THE INTERCHANNEL INTERFER- 
ENCE SPECTRUM 

Let Wa(f) be a finite power spectrum having limited total fluctuation 
for ^ / ^ » . Let the total power be finite so that the integral of w„(/) 
from / = to / = ^ converges absolutely. We define two auxiliary spec- 
tra 

M/), I/-/0I < e 

[ 0, l/-/o|>. 

w, = wM) - wXn (Al-l) 

and note that the autocorrelations corresponduig to these spectra must 
satisfy 

P^r) = PJr) - PM (Al-2) 

We consider the problem of transmitting the ensemble having the spec- 
trum Wb(f) through a system in which the input and output autocorrela- 
tion functions, ^i(t) and *2(r), respectively, are related by 

Here F(z) and its derivatives F'(z) and F"(z) are assumed to be finite 
continuous functions which exist over the range of 2 of interest. 

Let the output spectrum corresponding to the input spectrum Wb(f) 
be Wsif). We shall show that as « ^ the ^'alue of Wb(J) in the range 



INTERCHANNEL INTERFERENCE IN FM AND PM SYSTEMS 631 

\f~Io\ < t approaches 

4 ( {FIRJt)] - F'mR,(r)] COs27r/rrfT (Al-3) 

Ja 

When we multiply this expression by 2e and set/ = /o , we obtain the 
power appearing at the output of the system in an unloaded channel of 
width 2e centered on /o . Here we are not interested in values of Wsif) 
outside the range [ / ~ ,fo | < e. 

First we note several properties of autocorrelation functions. From 



R(t) = f wif) co8 2irJTdj 
Jo 



it follows that | R(t) \ ^ R(0). Also, if w(f) has limited total fluctuation 
in the interval (0, =0 ), the Riemann-Lesbesgue lemmas* and the absolute 
convergence of the integral for R(Q) show that Hir) = 0(1/t) as 7 — »■ » . 
Thus we may find positive numbers A, B, C such that for < r and 
any e less than some fixed value 

I lUr) I < A/r 

1 RAr) I ^ RM = { " wM) df < Be (Al-4) 

I RM I < C/t 

By the extended theorem of the mean the autocorrelation function 
corresponding to Wsif) is 

R,.(r) = F[R,(r)] = F[R„(t) - R^t)] 

= F[R„(t)] - F'[RM]R.(r) + r (Al-o) 

I r I ^ 2-'R:ir) I F"[R.{t) - eR^r)] [ < R:\r)D 

where ^ 6 g 1 and Z) is a positive number such that | F"{z) \ < D. 
Then 



Wb 



{/) = 4 f RM cos 2irfrdT = h- h-\- h 
Jo 

I^ = 4 r F[Rjr)] cos 2T/r rfr 
Jo 

(Al-6) 

I.,= 4 r F'iRJrMXr) cos 27r/r rfr 
Ja 



li = 4 [ r cos 27rfr dr 
Jo 



See, for example, Whittaker and Watson, Modern Analysis, 4th edition, p. 172. 
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Since F'[fl„(T)] ^ F'(0) + 5, where by the mean value theorem 

\s\ = \ R.{r) r'{dR.{T)] 1 < I K„(t) I D 

^ e g 1 

I1 may be written as the sum of two integralSj the second of which has an 
absolute value not greater than 

4 r \ R^{t)DRXt) Idr = 4D [ I R.(T)RXr) I dr 

Jo *'o 

+ 4D r I RAt)RXt) I dr < ADRMBeT 

+ 4Z) f ACr"' dr = 4D[RMBeT + AC/T] 

where T is an arbitrary number and we have used the inequalities 
(Al-4). Choosing T = e^^'^ shows that the last expression is 0(€^'^). 
Hence 

h = 4f'(0) r RXr) cos 2x/rdr + OC*'") 
Jq 

K(/), l/-/a| < e 

= 0(ei/^) + F'(0) 

[ 0, |/-/o|>. 

Therefore in the range \f — fo\ < e, which comprises the only fre- 
quencies of interest in the channel interference spectrum, 

h = 4F'(0) r Rjr) cos 27r/r dr + 0(e"') 
Jo 

By use of the inequalities f or | r | and | R^ir) \ we see that 

I /a I < 4 /" R,\T)DdT <W I BVdr -i- 4D r C'r-'dr 
jQ Ja Jt 

= 4D{B\^T + C/T) 

If we choose T = 1/e this expression is 0(e). 

When we collect our results and let e. become vanishingly small, we 
see that expression (Al-6) for Wb{J) approaches (Al-3) for frequencies 
in the range [ / — /o I < e. Thus in the limit as e — > we may use the 
autocorrelation function 

F[RM] ~ F'imM 
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to compute the interchannel interference spectrum. This is the result 

used in (1.22). 

Appendix II 

APPROXIMATE EVALUATION OF INTEGRALS OF A CERTAIN TYPE 

The problem of evaluating the integral G defined by (2.4) is quite a 
difficult one. Here we shall outline a method which often may be used to 
obtain an idea of the order of the magnitude of such an integral. 

Let F{u) be an even analytic function of u such that the major con- 
tribution to the value of 

I{a) = 2 [ F{u) cos audu = T F{u)e''''' du (A2-1) 

Jq J-„ 

comes from a saddle point on the positive imaginary w axis. Then the 
"method of steepest descents" suggests that an approximate value of 
/(o) may be obtained by the foUowing procedure. 

1. Betfiy) = F(iy) and plot z = d\[ogf(y)]/dy =^ f{y)/j{v) as a func- 
tion of y. 

2. Draw the horizontal line z = a. Suppose its first intersection with 
the curve obtained in step 1 is at y = yi , and let the slope of the curve, 
determined either graphically or by differentiation, be {dz/dy)y^ at ^i . 

3. Then 7(a) ^ [2ir/{dz/dy\f'j{y,)e-'"' (A2-2) 
It should be noted that (A2-2) cannot be used indiscriminately. Thus, 

it does not work well for F{u) = 1/(1 -\- u) because there is no saddle 
point on the imaginary w-axis. However, when it is applied to integrals 
of the type encountered in our study it appears to do fairly well, as 
Table 4.4 shows. 

Appendix III 

LEWIN's INTEGRAL 

Here we study the integral 

{e^ ainu _ ]^ _ ^y-1 gjj^ y^ ^pg audu (A3-1) 

00 

which occurs in several limiting cases in our work and which has been 
studied by Lewin^ for a = and a = 1. 

When the exponential term is expressed as a power series in (b sin u)/u 
and the result integrated termwise, we obtain 
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I{h,a) - t,Anbyn\ (A3-2) 

An = I ( ) COS au du 



2t 



2 (-)™C™"(n - 2m+ a)""' 



2"(n - 1)! ^ 

A„ -^ (67r/n)'^' exp [-3aV(2n)] 

where C^" = n\/m\{n-m)\ and the last term in the summation for An 
is the last one for which n - 2m + a is positive (assuming a 7^ integer) 
and for which m S n. When n ^ 2, 4 „ is a continuous function of a. 
Tables III A and III B were computed from (A3-2). 

When b is small, the first term in (A3-2) gives, for ^ a ^ 2, 

I(b, a) ^ &V(2 - a)/4 (A3-3) 

and when & is a large positive number the contribution of the exponential 
term in the region around )( = gives 

Kb, a) ^ (GT/bY" exp (b - ^) (A3-4) 

Lewin has given more carefid approximations for the a — and a = 1 
cases. 

When b is large and negative most of the contribution comes from 
around w = ±3Tr/2 where exp [{b sin u)/u] attains its largest values. 

It is found that 

I(-0, a) ^ 10.76 0-'" cos (4.49a)/''^'^-'-^°''''"' + R (A3-5) 

where = — & is a large positive number and K is a remainder term. The 
numbers in (A3-5) are related to the value Un = 4.493 . . . where (sin u)/u 
has a minimum. 

Computation shows that the value /3 = 8 is not large enough to make 
the leading term in (A3-5) a good approximation for /(— /3, a). In order 
to obtain a better approximation we write 

fw _ 

exp [— (3m ^ sin w] cos au du 

+ 2 / [— 1 + pu~^ sin id cos au du (A3-6) 

Jo 

+ / (/3m~^ sin uY cos au du 

Jy 
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-6, 



Table III A — e"7(6, o) for 6 > 



b 


c'' 


a - 


0.2S 


0.50 


0.75 


1.00 


1.25 





1 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.5 


1.649 


0.272 


0.241 


0.209 


0.176 


0.142 


0.107 


I.O 


2.718 


0.761 


0.685 


0.602 


0.511 


0.414 


0,314 


2.0 


7.389 


1.560 


1.440 


1.291 


1.117 


0.919 


0.713 


3.0 


20,08 


1.913 


1.801 


1.645 


1.448 


1.215 


0.968 


4.0 


54.60 


1.974 


1.888 


1.751 


1.566 


1.341 


1.098 


5.0 


148.4 


1.905 


1.844 


1.731 


1.571 


1.372 


1.153 


6.0 


403.4 


1,794 


1.751 


1.660 


1.525 


1.356 


1.166 


7.0 


1097. 


1.680 


1.649 


1.575 


1.463 


1.320 


1.157 


8.0 


2981. 


1.576 


1.552 


1.492 


1.398 


1.277 


1.138 



Table III B — I{b, a) for 6 < 



b 


n - 


.25 


.50 


.75 


[.0 


1.25 





0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


-0.5 


0.349 


0.300 


0.254 


0.210 


0.167 


0.125 


-1.0 


1,25 


1,06 


0.885 


0.723 


0.576 


0.432 


-2,0 


4,16 


3,41 


2.76 


2.20 


1.76 


1.34 


-3,0 


8,03 


6,37 


4.97 


3.88 


3.14 


2.46 


-4,0 


12.6 


i!,6{j 


7.23 


5.49 


4.55 


3.74 


-5.0 


17. K 


13.2 


9.40 


6.89 


5.93 


5.19 


-6.0 


23.6 


16.S 


11.4 


8.00 


7.25 


6.85 


-7.0 


30.0 


20.7 


13.1 


8.71 


8.48 


8.78 


-8.0 


37.2 


24. S 


14.5 


8.93 


9.59 


11.0 



where // is such that the quantity within the brackets in (A3-1), with 
h = — /3, is approximately (/3(r' sin w)V2 ^vhen H > y. For rough woi'k we 
may take /y = /3. The leading term in (A3-5) arises from the contrii:)ution 
of tiie region around n = Sir/'I to the value of the first integral in (A3-6). 
The contributions from the regions around u = 7ir/2, Uir/2, ... (if 
// is large enough) add to the \a\uc of R but they are generally small in 
comparison with the leading term in (A3-o). 

Thus we are led to approximate R in (A3-5) by the sum of the second 
and third integrals (expressed in terms of integral sines and cosines) in 
(AS-f)) with // = 0. When R is replaced by this sum, expression (A3-5) 
gives values for /> ^ — 8 which agree fairly well mth those in the table. 
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